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14, AfSTRACT 

During this second quarter, three activities were in progress: (1) the set-up 
of furnaces and photolithography related equipment at MSU, (J.D. Gassaway); (2) 
experiments on double layer metal at NASA/Huntsvllle (T.E. Wade); and (3) two 
dimensional modelling of MOSFETS (Q. Mahmood & J.D. Trotter). 

! The MSU T.C. Laboratory lias In place the furnaces, the photolithography facilitiet 

^ and metallization system to reproduce the NASA/Huntsville double layer metal process 

i the only missing equipment being the SIO^ vapox reactor. Two Pacific Western units ( 

I are in the process of being donated to MSU by Synertek. i 

* The double layer metal activity Initially emphasized wet chemistry techniques. 

Py incorporating the following techniques: (1) ultrasonic etching of the vlas, ! 

(2) premetal clean using a modified buffered HF, (3) phosphorus doped vapox, and (4) j 
extended sintering, yields of 98% were obtained using the standard test pattern. I 

The two dimensional modelling activity is In the middle of development. j 

Principle problems have stemmed from alternately, instability and too much computa trion j 
time to achieve convergence. Modifications are in progress, ^ 
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A. DOI'BUJ LAYKR MiriAL SINTERING EFFECTS 


Al. INTRODUCTION 

The increasing demands toward greater packing densities in LSI 
and VLSI make it imperative that multilevel metallization systems be 
daveloped. Although several reviews have been written the 
realization of reproducible and reliable results have yet to be 
forthcoming. The most common problem associated with double layer 
metal has been the inability to make electrical contacts between 
metal layers through via holes etched in the dielectric to provide 
electrical communication between metals. The problem is more acute 
when the number of vlas is large (>500) and are relatively small in 
size (<0.2 mil square). Another, less important problem is associ- 
ated with shorts between metal layers due to pinholes in the di- 
electric film, thin spots, or poor coverage of hillocks in the first 
level metal. An example of the use of double layer metal for cir- 

g 

cult application Is given in Figure 1. 

In order to study the effects associated with double layer 
metal structures having a large number of small vlas, a test pattern 
was generated consisting of a string of 560 vias. The via size was 
varied from 0.5 mil square to 0.2 mil square per string. Electrical 
measurements were made on the test pattern after initial sintering 
and subsequent h''at treatments in order to monitor contact resist- 
ance behavior and variations in yield. Metal interconnects between 
vlas were typically 0.7 mils wide although samples were also pre- 
pared and tested having interconnects of 0.4, 0.5, and 0.6 mil wide 
Al, with little or no variation in results. 
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COMPLETCO 
AL GATE C MOSFET 
DOUBLE LEVEL METALLIZATION 




Figure 1. An example of the use of double layer netal In the 
realization o1 an aluminum gate C-MOSFhF structure. 




All. PREPARATION OF TEST VEHICLE 


The steps Involved In preparing the test vehicle consist of the 

following. The starting material was 3-8 ohm-cm» (100) o' tented, 

n-type phosphorous doped silicon wafers. A field oxide was thermally 

grown In dry O 2 for 5 minutes, boiling H 2 O for 60 minutes and then 

dry 0, for an additional 15 minutes at 900*C resulting In an oxide 
K9 

thickness of 14 A . Prior to depositing the first layer metal, a 
cleaning step was performed, this consisted of a one-minute dip in 
dilute HF: II 2 O (1:10) at room temperature and a 10 minute rinse 
in dclonired water. The wafer was then dehydrated at 900"C in an 
N 2 ambient for 10 minutes. The first metal layer was d.c. sputtered 
structural grade A1 alloy 6061 of O.Spm thickness. The metal was 
patterned using conventional photolithographic techniques with 
Waycoat-31 negative photoresist. 

Prior to depositing the dielectric, scribe lines (0.7 mils wide) 
were etched around the test patterns through the thermally grown 
field oxide down to the silicon substrate. These lines were used 
as an etch end-point monitor in etching vias in the intermetal di- 
electric. 

The intermetal oxide was next deposited at 400“C using CVD of 

SiH^ (4% in Ar) and O 2 (and later ^2^5 approximately 3 mole per- 

K® 

cent) to a thickness of 8 A . Dielectric thickness was measured 
with an (laser) elllpsometer. Vlas were etched in this dielectric 
film with buffered HF using standard photolithographic masking. 

This etching was done first by dipping the wafers in a stirred 


solution of otchant» using the scribe lines as reference to deter> 
nine when etching was completely through the dielectric flln. In 
later proccsHon, this etching wan done ultr.iHonieally tn a totally 
enclosed container. 

The second layer metal was also deposited using dc sputtering 
of the same target used in the first layer metal and patterned. 
Typically a 15 minute aluminum sintering process at 470*C preceded 
the first testing of the completed test vehicle. 

AIII. EXPERIMENTAL RESULTS 

Approximately 50 wafers were processed. Each wafer had 200-220 
test patterns with via sizes ranging from 0.2 mils square to 0.5 
mils square. Table I represents a summary of results in terms of 
percent yield obtained for each via size as a function of fabrica- 
tion process. Process A through L are explained below. Among other 
variables it will be noted that each process is characterized by an 
etch time associated with etching the vias in the intermetal di- 
electric film. The scribe lines were used as reference. When all 
dielectric was etched from these scribe lines, this was defined as 
the 'break* time. Any additional etching beyond this break time Is 
measured in seconds. 

Process Definitions: 

A. Wafers were processed as described in section II. 

The width of the first level metal interconnect of 0.4 , 0.5, 0.6 and 
0.7 mils corresponded to test patterns having square vias of 0.2, 
0.3, 0.4 and 0.5 mils, respectively. The width of the second metal 


via Slza 


0.5 mils 


0.4 mils 


0.3 mils 


0.2 mils 


Ptfrccntngi' Yield 

Via Resistance for 560 Vlas > 20 Meg Ohm 
Process Identlficacltto 

A I B I C I I) i B I F I G I H I 1 i Jf K| L 


81 76 95 98 99 98 94 


94 92 98 98 


70 56 88 94 97 98 92 82 89 


97 97 


61 53 64 85 93 94 92 83 8/ 80 100 


40 18 64 82 91 94 93 89 72 67 98 98 


Table I- A 


Percentage Yield 

i^la Resistance for 560 Vlas < 10 Kll Ohms 
Process Identification 

Via Size j A I B I C 1 > I H 1 F | G | H j I | j| k! 


0.5 rails 17 59 86 97 |97 94 84 61 92 93 98 98 


0.4 mils 8 48 73 89 94 96 88 64 82 9o| 97J 97 


0.3 mils 0 30 45 64 89 92 1 88 74 1 85 80 100 100 

I I It 


0.2 mils I 0 14 13 26 44 59 64 52 38 


Table I-B 


98 

98 

97 

97 

too 

100 

98 

98 


Via Size 


0.5 mils 


0.4 mils 


0.3 mils 


0.2 mils 


Via 

Percentage Yield 
Resistance for 560 Vias < 1 

Kll 

Ohm 

A 1 

Process Identification 

b|c|dIe!f|gIh| 


j| k| l 


0 59 78 85 86 


0 46 47 66 72 72 68 j 54 «2 

I I I i 


0 28 85 36 47 55 58 56 85 80 


98 

98 

97 

97 


0 8 llO !15 


30 30 32 67 73 98 


Table I-C 


Table I Percentage yield as a function of processing with via 
size as a parameter for contact resistai.ee of 560 vias 
less than 20 Megohm, 10 Kllohm, and 1 Kllohro for tables 
1-A, I-B, and l-C, respectively. 







































































































































was 0*7 alia. Via atch time waa break plus 10 seconds. 

B. Wafers were processed as described In section IX. 

The natal interconnects on both levels wera 0.7 nils wide. (This is 
true for all proceHsen B through II). ' la etch ttmo was break plus 
10 seconds. 

C. Wafers processed as described in section II with via 
etch time in stirred dielectric etchant extended 20 seconds beyond 
break time. 

D. thru H. correspond to wafers processed as described in 
section II but with via etch time in stirred dielectric etchant ex- 
tended 40, 60, 80, 100 and 120 seconds beyond break, respectively. 
(Note that the values given for processes A through H represent an 
average of five different process runs of eight wafers for each run.) 

I. Wafers processed as described in section II but having 
a phosphorous doped intermetal dielectric of 3.1 mole percent (as 
determined by Auger spectroscopy) and with vias etched In stirred 
dielectric etchant for 20 seconds beyond break. 

J. Wafers processed as described In I above except via 
etching was done ultrasonlcally in a closed container for 20 seconds 
beyond break. 

K. Wafers processed as described In J above with the 
addition of a first level metal cleaning step through the vlas prior 
to depositing the second level metal. The metal cleaning solution 
consisted of an ethylene glycol-buffered HF-H 2 O solution for the 
removal of Al 20 ^ from the first level metal prior to deposition of 
second level metal^^. 


L. Wafers proccuHed ua described In K above but having 
undergone a one>half hour poet heat (aintering) treatment at 490*C. 

The reeulte of proeeaaee C through H are displayed graphically 
in Figure 2. The data are shovm in Tables lA, B» and C which give 
the percent yield for via contact resistances less than 20 meg-ohms, 
10 kil-ohms and 1 kll-ohm. It is desirable to have minimum contact 
resistance (i.e. less than 250 ohms) since as the number of vias 
increase from 600 to approximately 3000 for a VLSI circuit, the 
resistance will increase proportionately. Any pattern having a 
contact resistance greater than 20 meg-ohm Is assumed to be an 
open-circuit in Table I. 

Wafers processed as described in D in Table 1 were used for an 
Auger surface analysis of the first level metal through a 0.5 mil 
via Just prior to depositing the second level metal. The purpose of 
the Auger analysis was to determine what Insulating materials reside 
on the surface of the first level metal which prevented good low 
resistance ohmic contact between metal levels. The usults of this 
analysis indicated the following materials were present with an 
accuracy of ± 5%s 


A 1^03 

64% 

c 

31% 

N 

T.8% 

Mg 

1.3% 

Si 

0.6% 

F 

0.6% 

S 

0.7% 

0 (other than Al^O^) 

<0.1% 


For this contact, the amounts of Mg and Si were about the same as 


Yield (Percent) Yield (Percent) 


Yield vs. Etch Time 


0.2 Mil Vies 


< 20 M n- 


^ 10*^0 


< I’^fi 
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Etch Time Beyond Break (Seconds) 


Yield vs. Etch Time 
0.3 Mil Vlas 


M 

< 20"n 


< 10 *^ 5 ^ 


< 
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Etch Time Bevond Break (Seconds) 


Figure 2. Percent yield as a function of etch time beyond 
break with via size as a parameter. Intermetal 
dielectric etch accomplished with stirred B.O.E. 
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th« bulk A1 alloy so there did not seem to be any surface enrichment 
of Mg 2 Si which was suspected as a source of the high resistance. 

When about lOOA of this surface was removed by ion etching, the 
N,S,F, and C were greatly reduced, which suggested that they may be 
a source of interfacial contamination. 

The large percentage of carbon measured (31%) is believed to 
be a result of the photoresist leaving a carbon residue on the 
wafer. The composition of this photoresist could not be obtained 
from the manufacturer. 

It should be noted that the composition of the structural 
grade A1 6061 consist of the following impurities: 


Si 

0,4 - 0.8% 

Cu 

0.15 - 0.4% 

Fe 

0.7% 

Mg 

0.8 - 1.2% 

Mn 

0.15% 

Cr 

0.15 - 0.35% 

Zn 

0.25% 

Ti 

0.15% 


The primary reason for using this A1 alloy was for the prevention of 
hillock formation. 

Post heat treatmettt or sintering of a wafer exhibiting an 
initial poor yield (i.e. less than one kil-’ohm for a chain of 560 
vlas in series) can increase the effective yield by 500 to 700% as 
Illustrated in Figure 3 fur wafers lOH and 2W. These wafers were 
processed as in 'C above, having only a 15 minute 470°C sintering 
before initial testing. All addition sintering was done at 490^0 
in a nitrogen ambient. 






Average Contact Resistance for 560 Vias (Ohms) 
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The primary reason for this increase in yield with post heat 
treatment is due to the fact that the Al-to-Al contact resistance 
decreases with sintering) as demonstrated in Figure 4. Here) the 
average resistance for all test patterns having 0.5 mil square 
vias (approximately 55 of the total 212 chips) is plotted as a 
function of sintering time. As the time of sintering Increases, the 
number of test patterns having a contact resistance < 1 kll-ohm 
increases, hence the yield increases. 

On an individual chip basis, the contact resistance for 560 

vias per chip behaved as shown In Figure 5. For a given wafer 

whici) initially exhibits a relatively poor yield, the number of 

chips behaving in the manner shown for chips //IOC and lOW was 

approximately 51%, those behaving like chip //lOR was 27% and those 

like chip //lOT was 10%. The remaining chips either exhibited low 

initial contact resistance (■'300 obms) or were open-circuited 

(>20 meg-ohms). Those chips initially exhibiting a low contact 

resistance generally had their resistance lowered with increasing 

heat treatment, however a few chips showed a gradual increase in 

resistance with increased heat treatment time as shown in Figure 6. 

It is believed that this increase in resistance is due primarily to 

the consumption op A1 by SlO^ to lessen the A1 thlrkness. Floce 

at 490®C A1 reacts chemically with SIO 2 to form Al^O^ and free S1»A1 

o 15 

is consumed at a rate of 210A/hr. 

The most predominant treiul found In sintering of double level 
metal test patterns is domonstrnted In Figure 7, As Inferred 
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•% 



Time of Heat Treatment (Hours) 

Figure 5 . Contact resistance of 560 vias as a function of heat 
treatment time for several different chips of wafer 
//lOH at a temperature of 490“C. 




Contact Resistance of 360 Vias (Ohms) 






Cont 


Contact Realstance v. Time of H«at 
Treatment for Wafera #2 and tflO at 
A90*C. 



Figure 7 . Contact resistance trends as a function of heat treatment 
time for wafers H2 and //lO. 


from this plot of contact resistance as a function of post-lieac 
treatment for these two wafers, the average contact resistance 
decreases drastically for the first few hours of sintering. Because 
of this feature, the realisation of multiple level metal LSI having 
a large number of vlas (>2000) is a distinct posslbllty. 

One major problem involved in post heat treatment of double 
layer metal structures is that the dielectric (C.V.D. - Si02> has 
a tendency to crack or "craze" at elevated temperatures (>S00*C) . 

This cracking, as illustrated in Figure 8, is due to excess stress 

o 

on the Insulating layer . One method to lessen, if not-alleviate 
this problem, is to use a phosphorous doped dielectric of 3 to 5 mole 
percent. This increases the effective temperature for which sin- 

9 

tering can be done before cracking is observed . 

Using this phosphorous doped dielectric (also called phos- 
phosillcate glass or PSG) resulted in a much higher yield as shown 
in Table I and also a considerable less contact resistance as 
shown in Figure 9. Here the average contact resistance as a function 
of via size Is plotted for a test pattern with a doped dielectric 
layer for no heat treatment (curve B) and with 30 minutes post heat 
treatment at 490"C (curve A) . These results are compared to the 
best results obtained for the undoped case (curve C) , even after 
many hours of hunt treatment. The table in the figure indicates 
resistance and percent yield for each case presented. 
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A resistance laap of the wafer shown in ctrve B is given below: 

173 159 154 217 161 143 158 196 176 

170 171 ® 157 231 158 16 I 159 183 171 158 

166 163 168 158 166 226 162 142 162 171 1.67 161 164 

164 164 167 161 158 208 164 159 161 172 170 163 165 173 

168 170 170 160 159 212 163 161 163 228 169 164 • 172 164 

157 168 169 170 165 158 219 162 161 160 234 168 163 163 170 160 

164 170 170 171 164 156 213 161 163 159 221 173 159 163 170 165 

112 173 201 176 164 151 143 173 « 163 192 174 162 162 170 171 

164 175 375 181 167 159 405 172 162 158 181 167 161 164 167 165 

177 202 187 169 159 211 164 161 154 176 1''4 160 156 166 

176 216 184 171 161 228 169 158 154 » 162 156 154 165 

44 198 179 171 162 228 164 158 153 173 160 157 151 

260 181 171 163 231 166 159 154 172 160 153 

177 166 200 167 159 151 167 157 

Via size .4 .5 .2 .3 .4 .5 .2 .3 .4 .5 .2 .3 .4 .5 .2 .3 .4 


A IV. ANALYSIS AND CONCLUSIONS 

The following comments and conclusions may be derived from 
this investigation. 

1. ) The variation in width of the metal interconnects between 
vlas had little influence on yield, as expected. It may have had 

a slight effect on total resistance, especially for the 0.4 mil 
compared to 0. 7 mil wide interconnect, but this was not detectable. 
The major difficulty came in mask alignment of the smaller inter- 
connects. 

2. ) The extension of the via etch time of the dielectric film 


for both the dipped (stirred) etching and the ultrasonic etching 


Phot oml r I'OKraph (A) 



Phil t i>m i c’ roj;raph ( B) 

Figure 8. Pluitomicrogroph of test pattern at 300X n.agnlf icat ion 
Illustrating era.'Jng tendenry at t omperatt res >500®C 
and for fast pull from the ftirnace. Photomicrograph 
(A) shows crazing around a test p»id whereas (B) shows 
crazing in the test pattern itself. 
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Average Resistance of 560 Via Contacts (Ohas) 



Average Contact Restatancc V8» Via Sign 
For Wafer lA-1 Before and After Poet 


Curve 


iCurvc 



Heat Traatment 

• 

Via Size 

Avg. Resistance 

% Yield 

A .2 

206 

98 

.3 

170 

97 

.4 

160 

100 

.5 

150 

98 

B .2 

676 

98 

.3 

195 

97 

.4 

180 

100 

.5 

177 

98 

jC .2 

4.6M 

13 

1 :4 

289 

56 

241 

77 

\ .5 

237 

78 


Curve A 


'^urve C 





Curve B 


0.2 0.3 0.4 0.5 

Via Size (Mils) 

Average contact resistance vs. via size for wafer 1^14-1 
having phosphorous doped Intermetal dielectric before post 
heat treatment (curve B) and after 30 minutes at 490“C (curve A). 
For comparison with a wafer having an undoped dielectric, curve C 
represents such a wafer after undergone 8 hours of post heat 
treatment at 490®C. (Note contact resistance for 0.5 mil via 
of wafer in curve C prior to post heat treatment was 90 kllohms.) 


difficulty of roowving all of the diolectric film froa the via and 
tho formation of oxidca of the allcya ccntal.ied within the A1 all 
tend toward low yield, high reaiatance contacts. There are actually 
two approachea for aolving this problem, back sputtering and 
chemical etching. Back sputtering has the advantage that the oxide 
thickness can in principle be reduced to 2 cro the subsequent 
metallization is applied without exposure to air. The disadvant- 
ages include: possiblity of radiation Induced MOS damage; potential 

for contamination of tho sputtering system; redeposltlon of pre- 
viously sputtered materials especially from the substrate table; 
etc. The chemical etch procedure does not have these disadvantages, 
but the oxide thickness can never be reduced to zero.^^ 

The use of a 1:1:1 ethylene glycol: buffered -HF : H^O etch 

X IO 

of oxide without attacking the Al. (The 
second layer sputtered Al enn easily penetrate this oxide thickness 
to form good ohmic contacts.) This etch was done Just prior to 
second level metallization. 

6.) The application of post heat treatments can drastically 
increase yield by lowering total via contact resistance. There are 
several possible reactions which cou? occur to cause this temper- 
ature effect. Aluminum could react with aluminum oxide in the 
film between the layers to form various aluminum suboxides and 
thus break up the continuous layer of dielectric. Alternately, the 
recrystallization of aluminum in both of the levels could mechani- 
cally disrupt the thin native oxide layer, therefore forming alumi- 
num- to-aluminum contact. As yet, there is no evidence for this 


PiW£ BUMil HOT F&MCIT 


•peculatlcrii. 

Flguraa 4 and 5 Indicota that in some cases the contact 
resistance increases with post heat treatment time prior to the 
"diffusion" through this 'Interfacial contamination' by the two 
layers of aluminum. If this contamination were SIO^, then it 
could be sff J that there exists excess ionic silicon In the oxide 

12 

which becomes tied-up upon the application of post heat treatments^ . 
Since the contamlnailon is probably Al20^i one might say there exist 
excess ionized A1 in the oxide which reacts with the oxidizing 
species during the sintering process, thus decreasing the total 
ionic charge in the oxide layer causing the net resistance to 
Increase. 

7.) Testing the via test pattern can be tedious. For high 
resistance vlas the measured resistance was seen to be both light 
sensitive and current sensitive. As a result, all testing was 
accomplished in the dark. Ideally, very low testing potentials 
are applied in order not to break down a barrier layer thus turning 
a defective via into a good one. Also, theoretically, vlas should 
not be tested in series, since the entire applied voltage will be 
dropped across a defective via, probably causing it to break down 
and appear good. However, in our testing, all vlas were measured 
at once in series using a digital voltmeter. 

8). Recently, a new planar multilevel interconnection technol-- 
13 

ogy was introduced using polylroide films. Although this proce- 
dure appears to be very promising for future applications if taken 

at face value, it is felt there is still many facets of "magic" 
associated with It. 


B. OTHER DOUBLE LAYER METAL EXPERIMENTS 


Reviewing the literature has disclosed many investigations 
relative to the surface oxidation processes of aluminum, Including 
recent reported results by Bell Lab personnel which directly relates 
to double layer metal processing. These studies show that the 
aluminum surface is extremely reactive and grows AI 2 O 2 .nH 20 films 
rapidly in water based solutions to thicknesses which cannot be 
reduced during the sintering operation. Even the growth rate of 
on A1 in air at room temperature is reported to be in the 

O 

range of 30 A . This compares with the Bell Lab's estimate of 

O 

55 A of Al20^ which can marginally be reduced by sintering. 

Suspect processing steps are many; however, two stand out; 

(1) the via etch using buffered etch and (2) the initial oxidation 
during vapox deposition. 

Consequently, three sets of experiments have been started: 

1. Experiments aimed at measuring the AI 2 O 2 /AI interface using 
auger techniques. 

2. Experiments using plasma etch and/or nitride insulation. 

3. Experiments using in situ back sputtering. 

Only the last experiment has been completed. As was expected, the 
in situ back sputtering gave much improved results over the standard 
wet chemistry ~ the detail results being published elsewhere. 


C. TWO DIMENSIONAL MODELLING OF MOSFET 


As outlined In the last report key assumptions have been made 
in order to simplify the computation time: 

1. The mobile majority carriers are included in an infinites- 
imally thick layer of charge at the SI-SIO 2 interface. 

2. The current flow is then described by a one-dimensional 
equation, 

VT *'t 

‘ ■ f "n ‘ ^ ‘ " it, .01*1-. 

A C 

where 1 and s are channel current and charge per unit channel 
width. Eg is the tangential interface field and ;j^j and are the 
mobility factors accounting respectively for gate modulation and 
hot electron effects. 

3. The potential in the depletion region is defined by Pois- 
son's equation with appropriate fixed boundary conditions: a) fixed 

potential along the source drain and substrate, b) the horizontal 
component of the electric field under the source and drain and 

far from the channel is zero, and c) the potential in the channel 
is defined in conjunction with the current equation above plus 
the one dimensional Poisson equation in the gate oxide. 

Obvious limitations of this approach are associated with 
avalanche and bulk generated leakage currents - which normally are 
of secondary concern. Additionally, the mobile sheet charge assump- 
tion does not allow punch-thru or subthreshold current modelling; 
however, it is believed that this can be directly added to the 
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algorlchm at an appropriate time. The assumption of a one-dimen> 
sional electric field in the oxide is similar to the gradual channel 
approximation and, consequently, is of fundamental Importance. 

The initial algorithm utilized Iterative techniques in solving 
first Poisson's equation and then the current equation in the 
channel. The initial values of potential Inside of the boundaries 
were set at zero. Long computation times were encountered in 
solving Poisson's equation and instability occurred when Iterating 
through the current equation. 

The basic program has been modified to include the gate oxide 
region in the solution of the Poisson equation - in order not to be 
limited by the gradual channel approximation. It should be noted 
however, that routines based on the gradual channel approximation 
are being utilized in order to develop the Initial conditions for 
the main two-dimensional algorithm. In addition the program has 
been simplified by assuming constant channel mobility - a restric- 
tion which is to be removed at a later time. 

Emphasis is now being made on incorporating the following 
procedures: 1) Using the gradual channel approximation, calculate 

the channel node potentials. 2) Using the above with the given 
boundary conditions and the single dimensional d.c., solve for all the 
potentials in the oxide. 3) Using the alternating direction implicit 
method with the Newton-Raphson linearization technique, the Poisson 
equation is solved in the depletion region. 4) The polsson equation 
and the current continuity equation are simultaneously solved in 
the channel region. 5) Repeat 3) but extended through the channel 
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and oxide regions, again using the two-dimensional techniques. 

6) Repeat 4). 7) Repeat 5) and 6) until convergence is achieved. 
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